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Hadronic T decays offer the possibility of determining the strong coupUng a, at relatively low 
energy. Precisely for this reason, however, good control over the perturbative QCD corrections, 
the non-perturbative condensate contributions in the framework of the operator product expansion 
(OPE), as well as the corrections going beyond the OPE, the duality violations (DVs), is required. 
On the perturbative QCD side, the contour-improved versus fixed-order resummation of the series 
r*^ is still an issue, and will be discussed. Regarding the analysis, self-consistent fits to the data 

including all theory parameters have to be performed, and this is also explained in some detail. 
" ' The fit quantities are moment integrals of the T spectral function data in a certain energy window 

and care should be taken to have acceptable perturbative behaviour of those moments as well as 
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control over higher-dimensional operator corrections in the OPE. 
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1. Introduction 

The T lepton is the only known lepton heavy enough to also decay hadronically. In view of 
its mass of ^ l.SGeV, it provides an excellent laboratory for the investigation of low-energy 
QCD. On the other hand, the scale is so low that the inclusion of effects beyond perturbative QCD is 
required for precision studies. The non-perturbative effects appear in the framework of the operator 
product expansion (OPE) as vacuum-condensate terms, or contributions beyond the OPE, named 
duality violations, since they signal the breakdown of the quark-hadron duality picture of the OPE. 

The central experimental observable is the total T hadronic width 

rk" — ^ hadronsyrfr)! , , ^ , 

which was explored in the seminal theoretical analyses [2, 3, 4] as a means to determine the strong 
coupling tts . In subsequent years, the acquisition of large data sets by the ALEPH and OPAL col- 
laborations allowed the extraction of inclusive spectral T-decay distributions, the so-called spectral 
functions, and their separation into the light-quark (up and down) vector and axialvector channels, 
as well as the strange channel [5, 6, 7]. The additional experimental information can be employed 
to compute moments of the decay spectra which yield further data points to be compared to theory. 
In principle, also exclusive decay distributions are available and could be investigated, though they 
are less useful for basic QCD studies. 

On the theoretical side, recent years have seen several improvements in the description of 
hadronic T decays. The most important one was the impressive analytical computation of the per- 
turbative order correction [8], which has revived the interest in analyses from T decays. 
Generally, one might think that an additional order in the perturbative expansion would reduce 
the theoretical uncertainties. However, the theoretical expectation for the total hadronic width and 
moments of the decay spectra depends on the way that large logarithms appearing in the expan- 
sion are resummed by means of the renormalisation group (RG). The two most commonly used 
approaches ai^e fixed-order perturbation theory (FOPT), and contour-improved perturbation theory 
(CIPT) [9, 10], and the newly available ff{ccj) coiTection has made the dispersion in the theoretical 
prediction only more distinct. Hence, efforts in recent years have been devoted to better understand 
the origin of the differences in the RG improvement and attempts to single out the more reliable 
procedure [11, 12, 13, 14, 15, 16]. Those shall be reviewed in more detail below. 

The second main step forward has been a proper inclusion of violations of quark-hadron du- 
ality in the T sum rule analysis. The calculation of the theoretical side of the sum rule involves an 
integral over a circle in the complex energy plane with radius = or smaller (down to about 
1.5GeV^). In the region close to the physical axis with real, positive sq, the OPE breaks down and 
contributions beyond it become relevant. This should be clear, as the OPE is unable to directly de- 
scribe the hadronic resonance structure. Since a sound theoretical description of duality violations 
(DVs) is not available, a initial model was laid out in ref. [17], and later refined in refs. [18, 19]. 
A complete analysis requires the simultaneous, self-consistent determination of all occurring pa- 
rameters, that is, tti, the OPE condensate parameters, as well as the DV model pai^ameters. On the 
basis of the OPAL data, such an analysis was performed in refs. [20, 21], the discussion of which 
will be the second main topic of this writeup. 
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2. The anatomy of 7? T V 

Most suitable for a determination of are the T decay rates into light u and d quarks /?T,y /a 
via the vector or axialvector current, and the related moments, since in this case power corrections 
are especially suppressed. Theoretically, R^^ y can be expressed as 

2 



D>2 



(2.1) 



where Sew = 1-0201(3) [22, 23, 24] comprises electroweak corrections, 5^°' denotes the pertur- 
bative QCD correction, the 5,^^^^^ are quark mass and higher D-dimensional operator corrections 
which arise in the framework of the OPE, and 5yJ^ is the DV contribution beyond the OPE. Be- 
fore entering a more detailed discussion of the analysis, let us review our current knowledge 
regarding the anatomy of /?T,y+A> which takes the general form 



Rr.V+A = 3 iKdl'^EW [ 1 + + 



(2.2) 



The higher-dimensional OPE and the DV corrections have been lumped together into the non- 
perturbative correction 5^^^. To obtain the general picture, first we assume to be known. Evolv- 
ing the present PDG average as{Mz) = 0.1 184(7) to the T mass, yields Us^Mt:) = 0.3186(58). Even 
though the PDG average also includes oc^ from t's, removing this datum only has a small influence, 
so that it appears justified to employ the PDG value. 

To derive our expectation for d^^\ which gives the dominant correction to /?T,y+A> the theoret- 
ical expressions in FOPT and CIPT are briefly reminded. In FOPT the fixed renormalisation scale 
jj. = Mx is chosen, which results in [1 1] 



= £ «(M?)" £ kcn,kJk-i , (2.3) 



;(0) _ 

n=l k=\ 

where a{}JL^) = as{pi)/n, and c„^k are coefficients which appear in the perturbative expansion of 
the vector correlation function. At each perturbative order, the coefficients c„j can be considered 
independent, while all other Cn^k with k>2 are calculable from the RG equation. Explicit results 
for the analytically known c„.i can be found in ref. [8]. Furthermore, the 7/ are contour integrals in 
the complex 5-plane which for example are presented in [1 1]. 

In CIPT, on the other hand, an ^-dependent renormalisation scale is introduced, which partially 
resums higher-order effects, namely the ones due to the running of a^, leading to [9, 10] 



4? = Y^CnMM'i) (2.4) 



j(0) _ 

n=\ 

in terms of the contour integrals J"{MI) over the running coupling, defined as: 



J:{M^,) = 2^ j^-^{\-xf{\+x)a'\-M]x). (2.5) 



\x\ = \ 



In contrast to FOPT, for CIPT each order n only depends on the corresponding coefficient c„j. All 
contributions proportional to the coefficient c„.i, which in FOPT appear at all perturbative orders 
equal or greater n, are resummed into a single term. 
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Numerically, the two approaches lead to significant differences. Employing the value for 
as{MT) given above, one obtains 

5^2 = 0.2022 ±0.0069 ± 0.0030 = 0.2022(75) , (2.6) 
4? =0.1847 ±0.0048 ±0.0033 = 0.1847(58), (2.7) 

where the first error corresponds to the uncertainty in as{Mx) and the second to an estimate of 
higher orders through a variation of the coefficient cjj = 283 ± 283 [12]. Given the results for 
the value \Vud\ = 0.97425(22) [26], as well as /?t,v+a = 3.4671(82), which follows from 
BvA = 61.85(11)% and = 17.839(28)% [1], one can estimate 8^1^, with the finding 

S^P^PQ = -0.0086(80), 5^^4 = 0-0089(65). (2.8) 

This analysis shows that 5^^^ is expected to be < 1% and from this simple estimate compatible 
with zero at about la. Hence, the non-perturbative correction is much smaller than the perturbative 
one. Still, at the current level of precession, where the error on the PDG average for induces 
a shift in 5^°) of roughly 0.5%, the non-perturbative contribution becomes relevant. Before real 
progress in the determination of from T decays can be made, however, the difference between 
FOPT and GIFT has to be settled first, as it is certainly more important. 



3. Adler function at higher orders 

The question whether FOPT or GIPT provides a better approximation to 5'°^ hinges on the 
behaviour of the vector correlator, or equivalently the Adler function, at higher orders. To make 
progress in this direction, additional information beyond the analytically known orders has to be 
taken into account. An attempt towards this goal, based on a model for the Borel-transformed Adler 
function, was presented in ref. [12] and shall be described next. 

The perturbative expansion of the Adler function D{s) takes the form 

oo 

AT^Dis) = 1 +D{s) = £ c„,i a{sY . (3.1) 

n=0 

For the following it is slightly more convenient to utilise the function D{s) instead of D{s). Its 
Borel transform B[D\ [t) is defined by 

oo 

D{a) = J dte-'/"B[D]{t). (3.2) 



The integral D{a), if it exists, gives the Borel sum of the original divergent series. It was found that 
the Borel-transformed Adler function B[D] {t ) obtains infrared (IR) and ultraviolet (UV) renormalon 
poles at positive and negative integer values of the variable u = 9t/{4n), respectively [27, 28, 29]. 
(With the exception of m = 1, since there is no corresponding gauge-invariant D = 2 operator.) 

Guided by the large-j3o approximation [29], the influence of renormalon poles on the pertur- 
bative expansion should be as follows: asymptotically, that is for high orders, the perturbative 
expansion is dominated by the u = —\ UV pole that is closest to w = 0. At intermediate orders 
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some dominance of the low-lying IR poles (m = 2, m = 3) is observed, while very low orders are no 
longer dominated by a small subset of renormalon poles. This suggests the following ansatz [12] 

B[D]{u) = S[Dn(") +B[Df]{u) +B[Df]{u) +dl'' + d\°u, (3.3) 

which includes one UV renormalon at m = — 1, the two leading IR renormalons at m = 2 and m = 3, 
as well as polynomial terms for the two lowest perturbative orders. Explicit expressions for the UV 
and IR renormalon pole terms B[D^'^]{u) and B[D^^]{u) can be found in section 5 of ref. [12]. 

Apart from the residues t/^ ^ and d^^, the full structure of the renormalon pole terms is dictated 
by the OPE and the RG. Therefore, the model (3.3) depends on five parameters, the three residua 
d^^ , d^ and d^, as well as the two polynomial parameters and d\^. These pai^ameters can 
be fixed by matching to the perturbative expansion of D{s) up to ff{a^), whereby also the estimate 
for C5 J = 283 is required. The parameters of the model (3.3) then are found to be: 

d]''^ = -1.56-10-^, df=3.l6, df=-l3.5, (3.4) 

^/P0 = 0.781, t/f° = 7.66-10"^ 

The fact that the parameter d^'^ turns out small implies that the coefficient C2,i is already reasonably 
well described by the renormalon pole contribution, although it was not used to fix the residua. 
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Figure 1: Results for S^q (full circles) and 5^' (grey circles) at as(MT) = 0.3186, employing the model 
(3.3), as a function of the order n up to which the terms in the perturbative series have been summed. The 
straight line represents the result for the Borel sum of the series. 

The implications of the model (3.3) for 5^°^ in FOPT and CIPT is graphically represented in 
figure 1. The full circles denote the result for and the grey circles the one for , as a function 
of the order n up to which the perturbative series has been summed. The straight line corresponds 
to the principal value Borel sum of the series. The order at which the series have their smallest 
terms is indicated by the grey diamonds. As is evident from figure 1, FOPT displays the behaviour 
expected from an asymptotic series: the terms decrease up to a certain order around which the 
closest approach to the resummed result is found, and for even higher orders, the divergent large- 
order behaviour of the series sets in. For CIPT, on the other hand, the asymptotic behaviour sets in 
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earlier, and the series is never able to come close to the Borel sum. In the large-jSo approximation, 
this was already observed in ref. [30]. 

The superiority of FOPT over CIPT critically depends on the size of the u = 2 residue t/^. If 
this residue in real QCD would turn out substantially smaller than in eq. (3.4), CIPT could provide 
the better approach to the resummed series. Corresponding models have been studied in ref. [16] 
whose main aim was to investigate the perturbative behaviour of a large set of moments used in 
tty analyses from T decays within the Borel models. Nevertheless, there is no known mechanism 
in QCD to suppress the u = 2 residue and thus the behaviour that favours FOPT appears more 
likely. Other studies that implement information on the renormalon structure of QCD include 
refs. [13, 15], where also conformal mappings were applied to the series. 

4. Self-consistent fits, OPE and duality violations 

A self-consistent analysis of a, from t decays includes the simultaneous determination of all 
other parameters entering the game. These include QCD condensates as well as DV parameters. 
The required additional information can be obtained through the use of weighted integrals of the 
inclusive decay spectra up to an energy < M^, the so-called moments. Making use of the fact 
that the corresponding coiTclation functions are analytic in the complex i-plane except for a cut 
along the real axis, one can define the moments via [20, 21] 



where nj,^^ and n|/y^ are spin-1 and spin-0 mesonic correlators, and the particular case of eq. (2. 1) 
corresponds to /?^^^(M|) with the kinematic weight Wt:{s) = (1 -s/m1Y{\ + 1s/M]). 

Typically, in the past, moment analyses of T decay spectra were based on 5 moments to de- 
termine 4 parameters, a,, the gluon condensate (a^GG), as well as the D = 6 and D = 8 OPE 
corrections [5, 6, 7, 31]. These analyses, however, suffer from several deficiencies: the used mo- 
ments are all calculated at = M^. In ref. [32] it was then shown that the fit solutions are not 
stable under a variation of towards lower values. Next, because pinched-weights with zeros at 
s = so were used, it was assumed that DVs are negligible. Since the DVs enter differently in the 
different moments, this is potentially dangerous. And finally, the perturbative behaviour of some 
of the used moments is very bad [16], whence they should be avoided in the analysis. Besides, 
the ALEPH analyses of refs. [7,31] were based on data with an incomplete covaiiance matrix [33], 
and it remains open what is the impact of this shortcoming. 

The aim of refs. [20, 21] was to improve on all the above deficiencies. Firstly, only OPAL 
data [6] were used, employing the original set in ref. [20] and an updated set incorporating present 
day T branching fractions in ref. [21]. Next, DVs were included in the fit which required the 
use of one un-pinched moment particularly sensitive to DVs, the simplest choice being w{s) = 1. 
In multi-moment fits the additional moments were required to have good perturbative behaviour, 
which excluded moments sensitive to the gluon condensate, and to suppress D > 8 OPE terms. And 
finally, the 5o-dependence of the moments was included down to about ~ 1.5GeV^. Figure 2 




n: 



.(1+0) 
V/A 




(4.1) 
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Figure 2: Vector channel fit to the moments Ry^^ (^o)- J^sft panel: comparison of the experimental moments 
with the theoretical fit curves for FOPT in blue (solid) and CIPT in red (dashed). The (much flatter) black 
curves show the pure OPE parts without DVs. Right panel: comparison of the fit prediction for the vector 
spectral function with the experimental OPAL data. 

displays the outcome of the most basic fit with the single moment w{s) = 1 and just to the vector- 
channel spectrum. Also with regard to the OPE, this fit is cleanest as all OPE corrections are 
additionally suppressed. In ref. [21] the resulting was found to be: 

a,(MO = 0.325(18) (FOPT), a,(MO = 0.347(25) (CIPT). (4.2) 

The uncertainties in (4.2) turn out to be substantially larger than in previous a, analyses from the 
T. One reason is the use of OPAL data, but more importantly employing w{s) = 1. Naively one 
might think that the inclusion of further moments into the fit would reduce the error on a,. As 
demonstrated in refs. [20, 21], however, this is not so, because all moments are strongly correlated. 
Even combined fits to vector and axialvector with up to three weight functions only resulted in a 
mild reduction in the uncertainty on a^. 

In conclusion, it appears as if only better T spectral function data, either from a revised ALEPH 
analysis, or from the 5-factories BaBar and Belle, will help to improve the situation, together with 
a resolution of the CIPT versus FOPT controversy of treating higher perturbative orders. 
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